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Abstract 

We study the normal state spectral properties of the fermionic excitations in 
the Boson-Fermion model. The fermionic single particle excitations show a 
flattening of the dispersion as the Fermi vector k F is approached from below, 
forshadowing a Bogoliubov spectrum of a superconducting ground state. The 
width of the quasiparticle excitations near k F increases monotonically as the 
temperature is lowered. In the fermionic distribution function this tempera- 
ture dependence is manifest in a strong modification of n(k) in a small region 
below k F , but a nearly T independant n(k F ). 
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Intense theoretical work and controversy concerning the foundations for a novel non- 
Fermi liquid metallic ground state has followed the suggestion [|TJ that in high T c supercon- 
ductors the normal state might not behave as a standard Fermi liquid (FL). In attempting 
to fit the various anomalous thermodynamic and transport properties of these materials, 
the concept of a marginal FL was proposed |2| in which the imaginary part of the fermionic 
self-energy is supposed to vary like \u — e F \ rather than like (u — e F ) 2 as in an ordinary 
FL. Apart from this phenomenological approach to non-FL behavior, it is well established 
that ID systems show in fact a branch cut spectrum rather than poles for the single particle 
Green's function. For interacting Fermion systems such a non-FL behavior for ID is however 
not maintained if one goes to D > 1 [||]. 

Since possibly the non-FL properties of the high T c materials are a prerequisite for 
the large values of T c themselves Q and also since the concept of a non-FL in itself is of 
fundamental theoretical importance, it is of great interest to study models which exhibit 
transitions between FL and non-FL behavior as some parameter of the model is varied. 

In anticipation of a breakdown of FL behavior for the Anderson impurities lattice model 
with finite range interaction, the single impurity problem has been studied ||. When the 
repulsive interaction is increased beyond a certain threshold a change-over from FL to non- 
FL is indeed observed due to the coupling to the ensemble of the orbital channels rather 
than to only the one which hybridizes with the impurity orbital. 

There are so far few attempts to study models in which such effects are active in a 
cooperative fashion and for D > 1. Khodel and Shaginyan || have recently proposed that a 
3D system with long range Coulomb interaction shows a breakdown of FL behavior with the 
development of a horizontal inflexion point in the one particle excitation spectrum at the 
Fermi momentum k F . Such a feature leads to a restructuring of the Fermi filling and hence 
of the Fermi distribution function at T = resulting in a finite slope (rather than a jump) at 
k F . This phenomenon is very reminiscent of what happens in a BCS superconductor with a 
Bogoliubov like excitation spectrum. It has since been shown , that the Khodel- Shaginyan 
state is an artifact of the Hartree Fock analysis and that it is indeed unstable towards a 
superconducting ground state. Nevertheless it was recognized |7j that such inherent features 
like potentially flat portions of the Hartree Fock excitation spectrum should be very sensitive 
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when including particle - particle scattering. Quasi particle broadening should then become 
a dominant feature and possibly lead to results which resemble in nothing an ordinary FL. 

On the basis of a detailed analyses of the anomalous lattice properties of high T c materials 
we have recently attempted to argue for the applicability of the Boson-Fermion model to 
high Tc superconductivity |J. As we shall see in this letter, a detailed analysis of this 
model shows features which were recognized as potential ingredients for the breakdown of 
FL behavior |7] in connection with the study of the Khodel-Shaginyan model. The Boson- 
Fermion model was initially introduced by one of us |J in order to provide a heuristic 
description of intermediate coupling electron-phonon systems and to study the resulting 
superconducting properties. This model was then extensively studied (although in a slightly 
generalized form) in connection with high T c superconductors, in the superconducting state 
||10||.It is therefore compelling to examine its normal state properties and to compare them 
to the experimentally established anomalous properties - which are the real testing grounds 
for any theory of high T c superconductivity. 

The Boson-Fermion model describes systems with both localized Bosons with charge 2e 
and itinerant Fermions with charge e on a lattice. An exchange of a Boson with two Fermions 
with opposite spins is assumed and the conservation of charge requires a common chemical 
potential for the two species of charge carries. This model exhibits a superconducting ground 
state with a BCS like gap in the Fermion excitation spectrum and collective excitations of 
the Bosons. What is new in this model as compared to a standard BCS systems is, that well 
above T c a pseudo-gap in the density of states develops which ultimately is expected to open 
up into a true gap at the critical temperature T c . It is this regime of temperature above T c 
which is of interest to us here and where clear non-FL properties manifest themselves. 

In this letter we shall study the Boson-Fermion model 

H = (zt-fi^c+Ci, - t ]T ctc la + (A B -2fi)Y,bfh 

+ "Et^i^t + C 1 ) 

i 

for a square lattice, t denotes the electron hopping integral, z the number of nearest neigh- 
bors, A s the atomic level of the Bosons (bi-polarons) and v the strength of the Boson- 
Fermion exchange coupling. The commutation relations for the Fermion and Boson operator 

3 



are given by {c ia ,Cj a r} = 5ij5 aa i and b i ,b'j' = 5^ respectively. We evaluate the excitation 



spectrum of the Fermions in the normal state in a fully self consistent way [IT| to second 
order in the exchange interaction. The expressions for the self energies for Fermions and 
Bosons are: 

„,2 



E jF (k, ov, 



N 



G F (-k + q,^ m - LO n ) G B (q.,uJ„ 



E B (q,u; m ) = — Yl G F (-k + q, -u n + u m ) G F (k,u n ) 

k,w n 

This leads to the set of equations for the Fermion and Boson Green's functions 



(2) 



G F (k,u n ) = [iu n - e k - E F (k,w n )] 1 



G B (q.,u m ) = [iuj m - E - E s (q,u; m )] (3) 

which are determined selfconsistently, together with the expressions for the self-energies, 
Equ.(§). k and q denote the momenta, u n and oj m the Matsubara frequencies for Fermions 
and Bosons respectively and N the number of sites. 

The unperturbed Fermion dispersion including the chemical potential is given by ek = 
£k — A 4 ; £k — t(z — 5 e * k ^)) with 5 denoting the vectors linking nearest neighbor lattice 
site. The unperturbed Boson energies are given by Eq = A B —2fi; the factor two in front of 
the chemical potential taking account that each Boson is constituted of two Fermions. 

The self-consistent coupled equations ||) are solved by an iterative procedure in which 
G F (k, u n ) and G_B(q, u) m ) are evaluated for a set of Matsubara frequencies u n = 27ik B T(n+^) 
for —100 < n < +99 and u) m = 2iik B Tm for —100 < m < +100. As usual we compute 
the difference between the full and bare Green's functions, so that only a small number of 
Matsubara frequencies is necessary. We restrict ourselves in the present study to summing 
the k and q vectors over a two-dimensional Brillouin zone with a set of 41 x 41 equally spaced 
vectors for the Bosons as well as the Fermions. Convergency of the iterative solutions of these 
self-consistent equations is obtained relatively fast for temperatures down to T = 0.0085 in 
units of the Fermionic band width 8t. The solutions for the Fermion and Boson Green's 
functions in terms of the Matsubara frequencies are then analytically continued to the real 
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frequencies axis and into the lower half plane using a standard Pade approximants procedure 
in order to obtain the poles of the retarded Green's functions and hence the excitation spectra 
for the Fermions and Bosons. 

In Fig.(|l|) we plot the Fermion spectral function v4i?(k, u) = —21m Gp(k, iu n = uj + i0 + ) 
as a function of Fermion momentum k which clearly shows how the spectral weight is split 
and redistributed in the region near k F ; k F being defined by T> F (k F ,0) + ek F = 0. This 
behavior manifests itself in the Fermionic density of states Pf(w) in form of the appearance 
of a pseudo-gap which deepens with decreasing temperature (Fig.|TJ) . Throughout the present 
work we assume as parameters: A s = 0.4 and v = 0.1 in terms of the bandwidth, v is chosen 
in such way that for realistic temperatures (a few hundred degrees Kelvin) we can expect 
noticeable effects due to the Boson-Fermion exchange coupling. The total average number 
of particles per site, n tot = n F + 2n B = £ k(T < c t.,cr c k,o- > +jj E q < b+b^ >, Bosons and 
Fermions included, is taken to be unity. 

The spectral functions of the Fermionic single particle excitations near k F have not a 
simple Lorentzian behavior and are best described by the continued fraction Pade approx- 
imants for Gp(k,u). For practically all k vectors around k F , one obtains essentially three 
poles whose spectral weight adds up to essentially unity. These three poles are given by 
a cosin like dispersion - indicative of the unrenormalized Fermion excitation spectrum ek 
- and two flat dispersionless branches separated by an energy equivalent to fi. Solving 
the Hamiltonian, Equ.([IJ) in the atomic limit (t — *■ and A s < 0) we find that for states 
containing two charge carriers per site there are two distinct eigenstates having energy 



e± = A B /2 ± y (A s /2) 2 + v 2 — 2/x, e + ~ corresponds to a state where the charge carriers 
exist predominantely in form of on site bosons - a bonding state -, e_ ~ —2// ~ |A S | corre- 
sponds to a state where the charge carriers exist predominantely in form of two uncorrelated 
fermions - an antibonding state. 

The real (o>k) and imaginary (— 7k/2) part of these poles are presented in Fig.(Q.a, b) 
and the modulus of their residues is given in Fig.(^[c) as a function of Fermion momen- 
tum. The residues are in general complex numbers whose imaginary part is biggest for the 
modes with largest damping 7k- Complex residues in general indicate strong interference 
between poles, meaning that the single particle excitations are no longer described by a 
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simple pole of the Green's function as in a classical FL. In the description of deviations from 
FL properties the frequency dependence of the imaginary part of the Fermion self-energy, 
r(k, u) = — 2im£^(k, iu n = uo + i0 + ), is of capital importance. For a classical FL we 
have r(k, u) ~ aT 2 + b(oj — fi) 2 . In the Boson-Fermion model the deviations from such a 
FL behavior are signifiant as can be seen from Fig.(Q) where we plot r(k F ,u;) for various 
temperatures. As in a FL, r(k, to) turns out to be very weakly dependent on k. r(k F ,u;) 
at finite temperature has a minimum which occurs for slightly negative frequencies. The 
reasons for that are at present not understood. As the temperature decreases this mini- 
mum shifts towards uo = but r(kp, 0) ~ T~ a increases with decreasing temperature (with 
a = 0.086 for 2D and a = 0.60 for ID) which is exactly the opposite of what occurs in a 
standard FL. All these manifestly non-FL properties of the Boson-Fermion model show up 
in the Fermion distribution function n(k) only in a relatively minor fashion. n(k) turns out 
to be practically T independent at k F , (all lines of n(k, T) cross within a very small regime 
arround k F ) but increases substantially with decreasing temperature for k just below k F . 

The breakdown of FL properties characterized by the disappearance of well defined 
quasiparticles near the Fermi energy is found to go hand in hand with the appearance of well 
defined itinerant excitations of the intrinsically (bare) localized Bosons due to a precursor 
to superfluidity. Such a precursor induced coherence of intrinsically localized Bosons has 



recently been studied by us |12| on the basis of the ID Boson-Fermion model. 

As the temperature is decreased the low q vector excitations, manifest in the poles of 
Gb(<1, uj), acquire the spectrum lu^ = h 2 q 2 /2m B (T) with a temperature dependent mass 
which decreases upon decreasing the temperature. We estimate the decrease of m B {T) to 
be given by m B = 24.2, 25.7, 51.3 and 195.6 for T = 0.0085, 0.01, 0.02 and 0.05. 

The present analysis of the Boson-Fermion model shows that for ID as well as for 2D 
the FL properties are totally destroyed. This is due to strong superconducting fluctuations 
setting in well above T c as can be seen from the pair correlation function - essentially 
determined by the Bose Green's function Ge(q, uj) - which has perfect Lorentzian line 
shape well above T c . From the static part of Ge(q, oS) we can estimate the coherence 
length, determined by < &q& q >— a(l + ^q 2 )^ 1 which varies as £ ~ (T — T c )~ u (with 
v ~ 1.25 and T c ~ for 2D and v ~ 1.1 and T c ~ for ID. The precursor effects of the 
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onset of joint superfluidity of the Bosons and superconductivity of the Fermions can also be 
seen in how the Hugenholtz-Pines theorem (G^ 1 (q = 0, lo = 0) = at T c ) is approached i.e. 
E + i: B (0, 0) ~ (T-T c .)- a with a ~ 2.7 for 2D and a ~ 2.4 for ID. These precursor effects 
forshadow a Bogoliubov like spectrum in which flat portions of the Fermionic excitations 
spectrum appear near k F which however are heavily overdamped. As a consequence, a 
pseudo-gap in the Fermion density of states begins to open up well above T c . It gradually 
deepens and is expected to become a true gap at T c . At present our numerical analyses of 
the normal state does not permit to prove this but our previously mean field results as 
well as preliminary RPA results for the superconducting state suggest such a behavior. 

The fully selfconsistent scheme employed in the present study corresponds to the sim- 
plest approximation in such a treatment. We believe nevertheless that the results obtained 
here should be qualitatively correct since vertex corrections to the selfenergies are of or- 
der f 4 smaller than the selfenergies taken into account here. Moreover the Hugenholtz- 
Pines theorem in the limit T — > T c is fulfilled, satisfying the selfconsistent equation 
E - S B (q = 0, to = 0) = for T = T c = 0. 

At this points the Boson-Fermion model should be compared with the negative U < 
Hubbard model for which a two particle resonant state appears inside the two particle 
continuum. Contrary to the Boson-Fermion model, these Bosonic states are however heavily 



overdamped in the U < Hubbard model [TJ]] and thus are not expected to control the 
superconducting state via their Bose Einstein condensation. The absence of well defined 
bosonic excitations does not cause single Fermion excitations to become overdamped and 
the U < Hubbard model, in spite of many similarities with the Boson-Fermion model thus 
remains a standard FL in the intermediate coupling regime i.e. U ~ 4t |T4| . 

Concerning the applicability of the Boson-Fermion model to high T c superconductors, it 
shows at least one predominant feature which a series of experiments suggest; namely the 
existence of a pseudo gap and non-FL behavior in the normal state. NMR, Knight shift, 



susceptibility measurements, optical conductivity, specific heat, etc [15| all suggest such a 
pseudo-gap in the underdoped samples. The behavior of the width of the fermionic spectral 
function, measured by photo-emission, favors however more the picture of a marginal FL 
than a completely destroyed one as discussed in this letter. In this connection we should 
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however remember that the Boson-Fermion model examined here is the prototype of this 
model which does not take into account any details of the cristallin lattice and electronic 
structure of real materials. Further work on this matter, also including shortrange spin 
correlations in the Fermionic subsystem, will eventually elucidate these questions. 
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FIGURES 

FIG. 1. Density of states and the spectral functions of the electrons for a set of k vectors, 
k = k x = k y = |j- x [0...11]. The results for additional equidistant k vectors was obtained by an 
interpolation of the selfenergy. In the inset we show the deepening of the pseudo-gap, centered at 
energy 0, as the temperature decreases, for three characteristic temperatures. 

FIG. 2. a) The real part of the poles of the Fermion Green's function Gp(k, a;) for k vectors 
along the diagonal of the Brilloin zone k = k x = k y = |j x [0...20]. b) The imaginary part divided 
by the real part of the poles of Gf(k,o;). c) The modulus of the residues of the poles of Gp(k, uj) 

FIG. 3. The imaginary part of the electron self-energy for k x = k y = k F for various tempera- 
tures. r(k F ,0) increases with decreasing temperature. 
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